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This report discusses several topics in both top quark physics and QCD at an International Linear Collider (ILC). Issues 
such as measurements at the tt threshold, including both theoretical and machine requirements, and the determination 
of electroweak top quark couplings are reviewed. New results concerning the potential of a 500 GeV e + e — collider 
for measuring Wtb couplings and the top quark Yukawa coupling are presented. The status of higher order QCD 
corrections to jet production cross sections, heavy quark form factors, and longitudinal gauge boson scattering, needed 
for percent-level studies at the ILC, are reviewed. A new study of the measurement of the hadronic structure of the 
photon at a 77 collider is presented. The effects on top quark properties from several models of new physics, including 
composite models, Little Higgs theories, and CPT violation, are studied. 
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1. INTRODUCTION 

The precision study of both new and already discovered particles will be a major component of the experimental 
programs at both the Large Hadron Collider (LHC) and a future International Linear Collider (ILC). The past two 
decades in particle physics established the importance of this precision physics program. The outstanding success 
of the Z-pole program at LEP and SLC elevated the global fit to the precision electroweak data into the primary 
experimental constraint on extensions of the Standard Model (SM). When combined with input from the Tevatron, 
it probes energy scales far beyond the kinematic limits of current colliders. 
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The enormous production rates for top quarks at future colliders, reaching 10 7 tt pairs in a 10 fb" 1 year at the 
LHC and 10 5 in a 100 fb -1 year at a 500 GeV e + e~ collider, will allow a similar program studying the top quark 
to be pursued. Rare decays of the top quark, deviations from its chirality structure in the SM, and its electroweak 
couplings will be studied. A high precision measurement of the top quark mass will greatly reduce the uncertainties 
in important electroweak parameters. A variety of work is needed for this program to be successful, including the 
precision calculation of top quark properties in the SM, the determination of experimental capabilities for performing 
measurements, and finding the most likely deviations predicted by models of new physics. The calculation of higher- 
order QCD corrections to top quark cross sections can have an important effect on top quark physics; for example, 
the threshold corrections to the process e + e~ — > ttH increase its rate by a factor of two, drastically increasing the 
sensitivity of a 500 GeV ILC to the top quark Yukawa coupling (see Section 2.5 of this report). The precision needed 
for the tt threshold scan imposes strong requirements on the monitoring of the luminosity spectrum, which are 
discussed in Section 2.7. Detailed analyses can reveal surprising experimental possibilities, such as the measurement 
of the Wtb coupling below the tt threshold (see section 2.1 of this report). The study of the predictions coming 
from models of new physics show that the top quark can be a powerful discriminator between different theories; for 
example, composite theories such as the Randall-Sundrum model predict shifts in the coupling of right-handed top 
quarks to the Z, while Little Higgs models generically predict shifts in the left-handed top couplings (see Sections 4.1 
and 4.2 of this report). 

In addition to the study of the top quark, precision programs studying the Higgs boson or Higgs mechanism, the 
W and Z bosons, and other new particles discovered will be possible. To fully utilize the percent-level experimental 
precisions for each of these programs, higher order electroweak and QCD corrections in the SM must be included. 
For example, the two-loop QCD corrections to the g — 2 of the &-quark reach the 2 — 3% level, and must be included 
when studying b production during a Giga-Z phase at the ILC (see Section 2.3 of this report). Benchmark processes 
for studying a strongly-interacting Higgs sector such as VlVl — > receive QCD corrections reaching 10 — 20%, which 
can mask the effects of new physics if not taken into account (see Section 3.2 of this report). Finally, fundamental 
properties of QCD such as the running of the strong coupling constant and the hadronic structure of the photon can 
be studied with unprecedented precision (see Sections 3.1 and 3.3 of this report, respectively). 

In this report we discuss several issues in top quark physics and precision QCD, with a focus on the physics 
program at the ILC. Section 2 discusses precision studies of the top quark at the ILC. Important issues such as the 
measurement of the tt threshold cross section are reviewed, and new results such as the study of the Wtb coupling 
below the tt threshold and the measurement of the top quark Yukawa at \fs = 500 GeV are presented. A preliminary 
study on the precise determination of the average beam energy and luminosity spectrum at the ILC, required for the 
tt threshold scan, is presented in section 2.7. Section 3 discusses important QCD physics that can be performed at 
the ILC. Higher order QCD corrections needed for several important measurements are discussed, and new results 
for the determination of the hadronic structure of the photon are presented. Section 4 studies modifications of top 
quark properties in several models of new physics, and analyzes the potential of the ILC to measure these shifts and 
use them to discriminate between different extensions of the SM. 

2. PRECISION STUDIES OF THE TOP QUARK AT THE ILC 
2.1. Measuring the Wtb Coupling Below the tt Threshold 

P. Batra, T. Tait 

A crucial test of the top quark's electroweak interactions is the strength of the left-handed charged current inter- 
action W-t-b. In the Standard Model, unitarity of the CKM matrix implies that gwtb ~ gVtb ~ 9, but in extended 
models, including the simple extension by a fourth generation of fermions, this interaction can differ significantly 
from the SM expectation. Currently, it is known that the W-t-b vertex is sufficiently strong that the dominant top 
decay is t — > Wb, but even an imprecise measurement is lacking. Single top production at the Tevatron and LHC 
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Figure 1: Inclusive rates for e + e~ — ^ W + bW~b as a function of the center-of-mass energy for gwtb = gsM (black solid), 
gwtb = 2gsM (blue dashed), and g W tb = gsu/2 (red dotted). 

will help fill this gap in our knowledge, and is expected to lead to a measurement at the 10% level, dominated by 
systematics [1]. 

Unlike other top measurements, a direct test of the W-t-b coupling is challenging at a 500 GeV e + e~ collider. 
A scan over the tt threshold region is expected to yield precise measurements of many top parameters in the SM, 
including the top mass, width and Yukawa coupling (see [2] and this report for projections), while above-threshold 
measurements may constrain anomalous, non-SM Lorentz structures [3]. Nevertheless, only an indirect measurement 
of the left-handed W-t-b coupling is offered from the tt threshold region, by inferring its value from the SM relation 
and a precise value of the top width. If, for example, there is a small non-standard decay mode of top, it will alter 
the width and distort the inferred coupling. It would be more desirable to have a direct measurement of W-t-b, by 
making use of a process which is directly proportional to it. Close to the tt threshold, sensitivity to the coupling 
is quite weak, because the rate is essentially the tt production cross section times the branching ratios for t — > Wb. 
Since we expect that the BR is very close to one, it does not in fact depend strongly on the magnitude of the W-t-b 
interaction. Meanwhile, single-top production above threshold, which is sensitive to the W-t-b coupling, is swamped 
by the tt background unless a collision mode is present [4] . 

Just below the tt threshold, the reaction e + e~ — » W + bW~b still occurs, through a mixture of non-resonant 
Feynman diagrams as well as through off-shell top quarks. At center-of-mass energies far enough below 2mt but still 
above m*, the rate is dominated by contributions from the virtual tt diagrams in a kinematic configuration where 
one top is on-shell and the other is off-shell. The rate becomes very sensitive to the W-t-b interaction, by virtue of 
the off-shell top [5]. This is illustrated in Figure 1, which plots the cross section as a function of energy for several 
values of gwtb, assuming a 175 GeV top mass and a 115 GeV Higgs mass. All analysis was performed using the 
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MadEvent package [6] at tree level. The lines asymptote to the same value at both ends of the energy spectrum, 
as on-shell ti production dominates close to threshold and graphs not involving top dominate far below threshold. 
Both of these extremes are independent of the W-t-b coupling. Thus, energies in between these two extremes are 
suitable to measure gwtb- We avoid the region very close to 2m t (despite its large rate), because the details of the 
transition from off-shell to on-shell do depend sensitively on the top width, which could obscure gwtb if there are 
non-standard decay modes of the top. Instead, we focus on the energy y/s = 340 GeV, where good leverage on this 
coupling appears to be attained with small dependence on the width. We will explore the interplay between gwtb 
and T t below. 

Here we restrict ourselves to simple cuts to model the acceptance. To that end, we require the jets to have 
Pt > 10 GeV and rapidities \y\ < 2. We assume W bosons can be reconstructed with little background and for 
simplicity assume perfect 6-tagging efficiency and no mis-tags. We improve the purity by requiring that one of the 
b quarks and one of the reconstructed Ws reconstruct an invariant mass within mt ± 10 GeV, though we do not 
assume the charge of either the b or the W can be determined. The dominant background that is independent of 
the W-t-b coupling comes from diagrams with an intermediate Higgs, which can be eliminated by subtracting events 
with bb that have an invariant mass close to the Higgs mass, once the mass is known. However, we do not impose 
such a cut in order to retain the most statistics possible. 

The number of events will depend strongly on the top mass, the Higgs mass, the top width and gwtb- It is 
expected that the ILC will determine the top and Higgs masses to order 100 MeV or better, which is enough to 
render the uncertainty in the rates from these parameters much smaller than the expected statistical uncertainties. 
The remaining dependence on the width and gwtb allows us to determine a combination of both these quantities. 
To illustrate the results, we assume 100 fb _1 collected at y/s = 340 GeV. In Figure 2 we present the contours of 
constant event numbers in the plane of gwtb and T t which reproduce the expected SM event rate of ~ 1500 events. 
Also shown are the contours corresponding to \a and 2a deviations from such a measurement (assuming that the 
SM rate is observed and considering purely statistical uncertainties). The result is the expected bound one would 
obtain on gwtb and Tt, which can be combined with the Tt from the above-threshold scan to extract gwtb itself (or 
alternately one can go to lower energies where the sensitivity to T t is less, though at the price of the loss of some 
statistics). From Figure 2, we see that assuming the width is measured with an uncertainty of 100 MeV, gwtb can 
be measured to the 2% level, which would represent better than a factor of 5 improvement compared to the LHC, 
and a major improvement in our understanding of the W-t-b interaction. 

Many improvements on these rough estimates are possible. Certainly a more detailed and exhaustive study of the 
background would be interesting, as well as more sophisticated study of the signal, including higher order corrections 
and theory uncertainties, and detailed modelling of the W decays and the observability of the b quarks. In particular, 
higher order QCD and EW corrections to the signal will be essential to include in a realistic analysis in order to 
obtain the desired accuracy in gwtb, but are not likely to strongly change our conclusions as to how accurately one 
will be able to measure the coupling. Finally, since the cross section is strongly energy dependent, it could also be 
beneficial to consider the utility of a number of smaller data sets at several different energies below threshold. We 
leave such refinements for future work. 

2.2. Probing Electroweak Top Quark Couplings at the ILC and the LHC 

U. Baur 

Although the top quark was discovered almost ten years ago [7, 8], many of its properties are still only poorly 
known [9]. In particular, the couplings of the top quark to the electroweak (EW) gauge bosons have not yet been 
directly measured. Current data provide only weak constraints on the couplings of the top quark with the EW gauge 
bosons, except for the ttZ vector and axial vector couplings which are rather tightly but indirectly constrained by 
LEP data; and the right-handed tbW coupling, which is severely bound by the observed b — > s-f rate [10]. 

At an e + e~ linear collider with y/s = 500 GeV and an integrated luminosity of 100 — 200 fb _1 one can hope to 
measure the ttV (V = 7, Z) couplings in top pair production with a few-percent precision [11]. However, the process 
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Figure 2: Curve corresponding to the SM rate and its la and 2<r deviations in the plane of gwtb and IV Also overlaid is an 
expected measurement of Tt from the on-shell threshold scan with an uncertainty of 100 MeV. 

e + e~ — ► ^y* / Z — > ££ is sensitive to both ££7 and ttZ couplings and significant cancellations between the various 
couplings can occur. At hadron colliders, tt production is so dominated by the QCD processes qq — > g* — > ££ and 
<7<7 — > that a measurement of the ££7 and ££.Z couplings via qq — > j* /Z* — > ££ is hopeless. Instead, the ££l^ couplings 
can be measured in QCD ££7 production, radiative top quark decays in ££ events (££ — > r yW + W~bb), and QCD ££Z 
production [12]. ££7 production and radiative top quark decays are sensitive only to the ££7 couplings, whereas tiZ 
production gives information only on the structure of the ££if vertex. This obviates having to disentangle potential 
cancellations between the different couplings. In this section we briefly review the measurement of the ££l^ couplings 
at the LHC and compare the expected sensitivities with the bounds which one hopes to achieve at an e + e~ linear 
collider. 

The most general Lorentz-invariant vertex function describing the interaction of a neutral vector boson V with 
two top quarks can be written in terms of ten form factors [13], which are functions of the kinematic invariants. In 
the low energy limit, these correspond to couplings which multiply dimension- four or -five operators in an effective 
Lagrangian, and may be complex. If V is on-shell, or if V couples to effectively massless fermions, the number of 
independent form factors is reduced to eight. If, in addition, both top quarks are on-shell, the number is further 
reduced to four. In this case, the ££l^ vertex can be written in the form 

r«^ 2 ,g, g -) = -ze{ 7 , (^(fc 2 )+75^(fc 2 )) + g (« + *)" (^v(^)+75^ 2 ))} , (1) 

where e is the proton charge, mt is the top quark mass, q (q) is the outgoing top (anti-top) quark four-momentum, 
and k 2 = (q + q) 2 . The terms i^(0) and F^ A (0) in the low energy limit are the ££V^ vector and axial vector form 
factors. The coefficients F^ v (0) and ^^(O) are related to the magnetic and (CP- violating) electric dipole form 
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Figure 3: The differential cross sections as a function of the photon transverse momentum for -y£u e bbjj production at the LHC. 
Part a) shows the SM signal and the various contributions to the background. Part b) shows the SM signal and background, 
and the signal for various anomalous ttj couplings. 

factors. 

In tiV production, one of the top quarks coupling to V is off-shell. The most general vertex function relevant for 
tiV production thus contains additional couplings, not included in Eq. (1). These additional couplings are irrelevant 
in e + e~ — > ti, where both top quarks are on-shcll. 

In e + e~ — > ti one often uses the following parameterization for the ttV vertex: 

Tf> 2 , q, q) = le j 7 „ (FUe) + l5 FMk 2 )) + (F?v(k 2 ) + lA(k 2 )) } • (2) 

Using the Gordon decomposition, it is easy to show that Eqs. (1) and (2) are equivalent for on-shell top quarks and 
that the form factors F^ A and F^ A (i — 1,2) are related by 

F 1V — — ( F 1V + F 2v) ) F 2V = F 2V i FlA = —FlA> F 2A = - ^2A ■ ( 3 ) 

The most promising channel for measuring the ttj couplings at the LHC is pp — > jiu^bbjj which receives contribu- 
tions from tij production and ordinary ti production where one of the top quarks decays radiatively, t — > Wb-f. In 
order to reduce the background, it is advantageous to require that both &-quarks are tagged. We assume a combined 
efficiency of e 2 = 40% for tagging both &-quarks. 

The non-resonant pp — ► W{— * lv)-fbbjj background and the single-top backgrounds, (tb-f + ibj) + X, can be 
suppressed by imposing invariant and transverse mass cuts which require that the event is consistent either with tij 
production, or with ti production with radiative top decay [12]. Imposing a large separation cut of AR(-f,b) > 1 
reduces photon radiation from the b quarks. Photon emission from W decay products can essentially be eliminated 
by requiring that m(jjj) > 90 GeV and mT{^]i> T ) > 90 GeV, where m(jjj) is the invariant mass of the jj*y system, 
and mT{(-^;i> T ) is the i"ii> T cluster transverse mass, which peaks sharply at raw . After imposing the cuts described 
above, the irreducible backgrounds are one to two orders of magnitude smaller than the signal. 

The potentially most dangerous reducible background is tij production where one of the jets in the final state fakes 
a photon. In Fig. 3a we show the photon transverse momentum distributions of the ti"f signal and the backgrounds 
discussed above. The tij background is seen to be a factor 2 to 3 smaller than the tij signal for the jet-photon 
misidentification probability (Pj^ 7 = 1/1600 [14]) used. 

The photon transverse momentum distributions in the SM and for various anomalous ttj couplings, together with 
the pt(j) distribution of the background, are shown in Fig. 3b. Only one coupling at a time is allowed to deviate 
from its SM prediction. 
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Figure 4: a) The differential cross sections at the LHC as a function of pr(Z) for (! I' ivbbjj final states. Shown are the 
SM predictions for tiZ production, for several non-standard ttZ couplings, and for various backgrounds. Only one coupling 
at a time is allowed to deviate from its SM value, b) The differential cross sections as a function of the missing transverse 
momentum for ^ T 6fe+4j production at the LHC. Shown are the SM predictions for tiZ production and for various backgrounds. 



The process pp — > tiZ leads to either I' l'~ Ivbbjj or £' £' bb + Aj final states if the Z-boson decays leptonically 
and one or both of the W bosons decay hadronically. If the Z boson decays into neutrinos and both W bosons decay 
hadronically, the final state consists of j> T bb + Aj. Since there is essentially no phase space for t — > WZb decays 
(BR(t — > WZb) ps 3 • 10~ 6 [15]), these final states arise only from tiZ production. 

In order to identify leptons, b quarks, light jets and the missing transverse momentum in dilepton and trilepton 
events, the same cuts as for ifcy production are imposed. One also requires that there is a same-flavor, opposite-sign 
lepton pair with invariant mass near the Z resonance, mz — 10 GeV < m(££) < mz + 10 GeV. 

The main backgrounds contributing to the trilepton final state are singly-resonant (tbZ + ibZ) +X (tbZjj, ibZjj, 
tbZlv and ibZiv) and non-resonant WZbbjj production. In the dilepton case, the main background arises from 
Zbb + Aj production. To adequately suppress it, one additionally requires that events have at least one combination 
of jets and b quarks which is consistent with the bb + Aj system originating from a ti system. Once these cuts have 
been imposed, the Zbb + Aj background is important only for pt(Z) < 100 GeV. 

The Z boson transverse momentum distribution for the trilepton final state is shown in Fig. 4a for the SM signal 
and backgrounds, as well as for the signal with several non-standard ttZ couplings. Only one coupling at a time 
is allowed to deviate from its SM prediction. The backgrounds are each more than one order of magnitude smaller 
than the SM signal. Note that varying Ff VA leads mostly to a cross section normalization change, hardly affecting 
the shape of the pt(Z) distribution. 

For the j> T bb + Aj [16] final state at least 3 jets with pt > 50 GeV and jfi T > 5 GeV 1 / 2 y/J^PT are required. The 
largest backgrounds for this final state come from ti and bb + Aj production where one or several jets are badly 
mismeasured, from pp — ► ttjj with tt — > t^vghbjj and the charged lepton being missed, and from tij production, 
where one top decays hadronically, t —> Wb — > bjj, and the other via t — > Wb — * rv T b with the r-lepton decaying 
hadronically, r — > hv T . 

In Fig. 4b we show the missing transverse momentum distributions of the SM ttZ — * p'rpbb + Aj signal (solid curve) 
and various backgrounds. The most important backgrounds are ttjj and tij production. However, the missing 
transverse momentum distribution from these processes falls considerably faster than that of the signal, and for 
■p T > 300 GeV, the SM signal dominates. 

The shape and normalization changes of the photon or Z-boson transverse momentum distribution can be used to 
derive quantitative sensitivity bounds on the anomalous ttj and ttZ couplings. For tiZ production with Z — * l' + £' , 
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Table I: Sensitivities achievable at 68.3% CL for the anomalous ttV (V = 7, Z) couplings FYv,a an d F?v,a °f Eq. (2) at the 
LHC for integrated luminosities of 300 ft -1 , and the ILC with yfs = 500 GeV (taken from Ref. [11]). Only one coupling at a 
time is allowed to deviate from its SM value. 



coupling LHC, 300 fb" 



[11] 



coupling LHC, 300 fb" 



[11] 



+0.043 
-0.041 
+0.051 
-0.048 
+0.038 
-0.035 
+0.16 
-0.17 



+0.047 
-0.047 
+0.011 
-0.011 
+0.038 
-0.038 
+0.014 
-0.014 



, 200 fb" 1 


AFfv 


, 100 ftr 1 


AF? A 


, 200 ftr 1 


AF 2 Z V 


, 100 ftr 1 


AF 2 Z A 



+0.24 
-0.62 
+0.052 
-0.060 
+0.27 
-0.19 
+0.28 
-0.27 



+ - 012 200 fb" 1 
-0.012 ' zw 10 



+0.013 100fb -i 
-0.013 ' 1UU 10 
+0.009 200 fb -i 
-0.009 ' zw m 
+0.052 10Q - -i 
-0.052 ' 1UU tD 



the A$ {£'£') distribution provides additional information [12]. In the following we assume a normalization uncertainty 
of the SM cross section of AW = 30%. 

Even for a modest integrated luminosity of 30 fb -1 , it will be possible to measure the ttj vector and axial vector 
couplings, and the dipole form factors, with a precision of typically 20% and 35%, respectively. For 300 fb -1 , the 
limits improve to 4 — 7% for F^ v A and to about 20% for F^ v A . At the SLHC, assuming an integrated luminosity 
of 3000 fb -1 , one can hope to achieve a 2 — 3% measurement of the vector and axial vector couplings, and a 10% 
measurement of F^y A , provided that particle identification efficiencies are not substantially smaller, and the reducible 
backgrounds not much larger, than what we have assumed. 

To extract bounds on the ttZ couplings, we perform a simultaneous fit to the pt(Z) and the A<&{£'£') distributions 
for the trilcpton and dilcpton final states, and to the f T distribution for the j> T bb + 4j final state. We calculate 
sensitivity bounds for 300 fb -1 and 3000 fb -1 at the LHC; for 30 fb -1 the number of events expected is too small 
to yield meaningful results. For an integrated luminosity of 300 fb -1 , it will be possible to measure the ttZ axial 
vector coupling with a precision of 10 — 12%, and F% V A with a precision of 40%. At the SLHC, these bounds can be 
improved by factors of about 1.6 (F% v A ) and 3 {Ff A ). The bounds which can be achieved for F^ v are much weaker 
than those projected for F^ A . As mentioned in Sec. 4, the Pt(Z) distributions for the SM and for F^ VA = —F^ V S A 
are almost degenerate. This is also the case for the A$> (£'£') distribution. In a fit to these two distributions, therefore, 
an area centered at AFf VA = -2F^ V S A remains which cannot be excluded, even at the SLHC. For Ff v , the two 
regions merge, resulting in rather poor limits. 

It is instructive to compare the bounds for anomalous ttV couplings achievable at the LHC with those projected 
for the ILC. The most complete study of ti production at the ILC for general ttV (V = 7, Z) couplings so far is that 
of Ref. [11]. It uses the parameterization of Eq. (2) for the ttV vertex function. In order to compare the bounds of 
Ref. [11] with those anticipated at the LHC, the limits on F^ V A and F% V A have to be converted into bounds on F^ V A 
and F^y A . Table I compares the bounds we obtain for F^ v A and F?f v A with those reported for the ILC in Ref. [11]. 
Note that only one coupling at a time is allowed to deviate from its SM value [11]. We show LHC limits only for 
an integrated luminosity of 300 fb -1 . For the SLHC, with 3000 fb -1 , we obtain bounds which are a factor 1.3 — 3 
more stringent than those shown in Table I. Thus, even if the SLHC operates first, and the $ T bb + 4j final state 
is taken into account, a linear collider will still be able to significantly improve the ttZ anomalous coupling limits, 
with the possible exception of Ff A . The ILC will also be able to considerably strengthen the bounds on F^ A and 
F^ A . It should be noted, however, that this picture could change once correlations between different non-standard 
ttZ couplings, and between ttj and ttZ couplings, are taken into account. Unfortunately, so far no realistic studies 
for e + e~ — > ti which include these correlations have been performed 1 . 

The LHC will be able to perform first tests of the ttV couplings. Already with an integrated luminosity of 



1 However, see Ref. [17] for limits on the CP-violating couplings. 
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30 fb -1 , one can probe the tt-f couplings with a precision of about 10 — 35% per experiment. With higher integrated 
luminosities one will be able to reach the few percent region. With the exception of Ff A , the ttZ couplings can only 
be measured with a precision of 15 — 50%, even at the SLHC. The ILC will be able to further improve our knowledge 
of the ttV couplings, in particular in the ttZ case. 

2.3. Two-loop Corrections to Heavy Quark Form Factors 

T. Gehrmann 

The international linear collider will produce large numbers of top-antitop quark pairs, thus allowing for precision 
studies of the properties of the top quark. These experimental precision studies require equally precise theoretical 
predictions, i.e. higher order corrections in perturbation theory. Up to now, the theoretical effort was focused on a 
precise description of top quark production at threshold (see [18] for a review), where QCD corrections to next-to- 
next-to-leading order (NNLO) in perturbation theory are known, while observables other than the total production 
cross section in the continuum are known only to next-to-leading order (NLO) accuracy. In this section, we present 
results on the virtual two-loop corrections to vertex functions involving heavy quarks, which are an important 
ingredient to the NNLO corrections to top quark observables in the continuum. 

The vertex function coupling an on-shell heavy quark antiquark pair to an external current can be decomposed 
into so-called form factors, whose coefficients follow from Lorentz invariance and symmetry properties of the current. 
For the vector and axial vector current (electroweak gauge boson), the vertex function contains, within QCD, four 
form factors (Fi, 2 , Gi. 2 ): 

= H) (vqF^s,™ 2 )^ + v Q ^-F 2 {s,m 2 )io-^( Pl + p 2 ) u 

+a Q G 1 {s,m 2 )Ylb + aQ-^G 2 (s,m 2 )j 5 (p 1 + p 2 Y 

The coupling of heavy quarks to Higgs bosons of positive and negative parity contains the scalar and pseudoscalar 
form factors: 

pi 

».~<^ = ~»™ [S Q F s {s,m 2 ) +iP Q F P (s,m 2 ) lb \. 

Here s is the invariant momentum squared of the external current and m the heavy quark mass. 

The two-loop corrections to the form factors are obtained by applying appropriate projections to the Feynman 
diagrams displayed in Figure 5. As a result, the form factors are expressed in terms of hundreds of different scalar 
integrals. These integrals are reduced to a small set of master integrals by means of the so-called Laporta algorithm 
[19] with the help of integration-by-parts identities [20] and Lorentz-invariancc identities [21]. The master integrals 
themselves were evaluated with the method of differential equations [21-24] in [25, 26]. The master integrals, and thus 
the form factors are represented as series in the regularization parameter e and expressed in terms of 1-dimensional 
harmonic polylogarithms up to weight 4 [27, 28]. 

We obtained the complete two-loop corrections for the renormalized vector [29], axial- vector [30, 31], scalar and 
pseudoscalar [32] form factors. Agreement was found with earlier partial results and with expansions around special 
kinematical points. 

An immediate physical application of the QCD corrections to the heavy quark form factors for arbitrary momentum 
transfer are predictions for form factors at zero recoil, the so-called static form factors. 

Measurements of these static form factors for heavy quarks could allow indirect constraints on new physics sce- 
narios. Recently, there has been considerable effort to determine the feasibility of such experimental measurements. 
Specifically, the couplings to photons and Z bosons have been studied in detail - both for heavy quark produc- 
tion at hadron colliders [12, 16, 33] and at a future high-luminosity high-energetic linear electron-positron collider 
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[11, 16, 34, 35]. At this workshop, this issue was revisited in great depth especially in view of discriminating new 
physics scenarios (see sections 4.1 and 4.2). 

The most prominent static form factor is the electromagnetic spin-flipping form factor: the anomalous magnetic 
moment, which is finite in the zero-recoil limit, and can be obtained from the results of the previous section [37]. It 
reads 

**<-o)-£ c, + 

with 

Fg$ = + 2 C2 (5 - 6 ln(2)) + 3 Cs) + C F C A + 3 ( 2 (-1 + 2 ln(2)) - | ( 3 

+C F T F (^--8 C 2 ) - y C F T F N, + C F f3 H^/ml), 
from which we can derive the static magnetic and weak magnetic form factor of a quark Q. We consider 

(Sr) 7 Q ^ f 2:.q (°) - v q z F ^ (°) • 

which correspond to the anomalous magnetic (MDM) and weak magnetic (WMDM) moments of Q. (Notice that in 
the literature the WMDM is often associated with F%q{s = m 2 z ).) Numerical values for t and b quarks are given in 
Table II. 

For the b quark an upper bound on its magnetic moment was derived in [38] from an analysis of LEP1 data, 
which, in our convention, reads \6(g — 2)^/2| < 1.5 x 10~ 2 (68 % C.L.). Comparing it with Table II we see that 
the QCD-induced contributions to the b quark magnetic moment saturate this bound, which implies that there is 
limited room for new physics contributions having the same sign as the QCD contributions to this quantity. At a 
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Table II: One- and two-loop QCD contributions, and their sums, to the anomalous magnetic and weak magnetic moments of 
the top and bottom quark, for different values of the renormalization scale fi. 

future linear collider [11, 35], when operated at the Z resonance, the sensitivity to this variable could be improved 
substantially, either by global fits or by analyzing appropriate angular distributions in bb and 667 events. 

As to the static form factors of the top quark, no such tight constraints exist so far on possible contributions 
from new interactions (see Ref. [12] for a review). These quantities are particularly sensitive to the dynamics of 
electroweak symmetry breaking. For instance, in various models with a strongly coupled symmetry breaking sector 
one may expect contributions from this sector to the static t quark form factors at the 5 - 10% level [118]. The 
QCD-induced anomalous magnetic moment and the QCD corrections to the axial charge of the top quark are of the 
same order of magnitude. Future colliders have the potential to reach this level of sensitivity. 

The form factors presented here have a number of applications, which will be addressed in future work. The vector 
and axial vector form factors contribute to the NNLO corrections to the forward-backward asymmetry for heavy 
quarks; they can also be used to compute the differential top quark pair production cross section at the ILC in the 
continuum, where top quark mass effects are still non-negligible. The scalar and pseudoscalar form factors enter 
the NNLO corrections to the decay of a Higgs boson into heavy quarks. These would become especially important 
for a very heavy Higgs boson, decaying into top quark pairs, where the form factors could be used for a differential 
description of the decay final state. 

2.4. Electroweak Effects at the ti Threshold 

A. Hoang 

The proper treatment of electroweak effects plays an important role in high precision measurements at the ti 
threshold. In the calculation of the total cross section, one can categorize electroweak effects into three classes. 

(a) "Hard" electroweak effects: This class includes electroweak effects related to the ti production mechanism itself 
or factorizable corrections to various matching conditions of the effective theory that is used to describe the 
nonrelativistic QCD dynamics of the top pair. In general these corrections are modifications to the hard QCD 
matching conditions of the effective field theory operators. They can be determined by standard methods and 
are real numbers. 

(b) Electromagnetic effects: Electromagnetic effects are relevant for the luminosity spectrum of an e + e~ initial 
state (beam energy spread, beamstrahlung, initial state radiation) and the tb final state (modification to the 
Coulomb attraction, contributions to hard electroweak effects). 

(c) Effects related to the finite top quark lifetime: Apart from the top decay (into Wb for the Standard Model) 
this class also includes interference contributions with processes having the same final state (W + W~bb) but 
only one or even no top quark at intermediate stages. This class also accounts for interactions involving the 
top decay products (" non-factorizable" effects). 

For processes involving highly energetic top quarks, where the momentum transfer typically involves scales much 
larger than T t , the class (c) effects are small corrections. On the other hand, for the top threshold dynamics all three 



PLEN0043 



2005 ALCPG & ILC Workshops - U.S.A. 



classes have to be accounted for at leading order, because the typical kinetic energy of the nonrelativistic top quarks 
is of the same order as the top quark width, 

E kin ~ m t v 2 ~ m t a 2 s ~ T t ~ m t a cm . (4) 

The relation in Eq. (4) gives the power-counting that has to be employed to systematically account for electroweak 
effects within the nonrelativistic expansion: 

v 2 ~ a 2 ~ a om . (5) 

Several analyses concerning electroweak effects belonging to the various classes have been carried out in the past. 
However, no coherent and comprehensive treatment that systematically accounts for all these electroweak effects 
beyond leading order currently exists. In the following a brief status report is given concerning top threshold 
production in e + e~ annihilation. 

The leading order (LL) electroweak effects belonging to class (a) describe the production mechanisms of the 
nonrelativistic top quark pair in e + e~ annihilation in the various possible spin and angular momentum states. Due 
to the power-counting, the one-loop hard electroweak corrections already contribute at next-to-next-to-leading order 
(NNLL) and involve the standard (real parts of the) one-loop electroweak corrections to the e+e~ — > process [39, 40]. 
There are also hard electroweak corrections that modify the QCD potentials [41], but due to gauge cancellations 
these corrections do not contribute at NNLL order. The dominant corrections caused by an exchanged Higgs boson 
through a Yukawa potential can also be considered as a class (a) contribution. 

The most important contribution belonging to the electromagnetic effects in class (b) is the luminosity spectrum 
that affects the cm. energy available in the e + e~ collision. While bcamstrahlung and the beam energy spread are 
machine-dependent and will have to be measured experimentally, the initial state radiation component is calculable. 
The luminosity spectrum leads to sizable smearing of the cross section and contributes at LL order. Electromagnetic 
effects also modify the QCD potential through photon exchange between the top pair and can lead to hard electroweak 
corrections belonging to class (a). In present analyses QED effects arc only accounted for through the luminosity 
spectrum. In particular there is no coherent treatment of the calculable QED effects that systematically accounts 
for the effects of initial state radiation, the Coulomb corrections and the hard QED effects. 

The leading order effect belonging to class (c) is the top decay width, which makes the perturbative treatment 
to the strong tt dynamics at all possible [42, 43]. Within the nonrelativistic effective theory used to describe QCD 
effects the top quark width can be implemented through an imaginary mass shift in the top quark propagator 



where 5m t is related to the top quark mass definition that is used. This effect can be incorporated into the forward 
e + e~ — ► e+e~ scattering amplitude results for stable top quarks by simply shifting the cm. energy into the complex 
plane: yfs — > y/s + iT t . Using the optical theorem one can obtain the total cross section. Beyond leading order the 
electroweak effects belonging to class (c) can be systematically incorporated into the effective theory by accounting for 
absorptive parts related to the top decay final states in the matching conditions [44] of the effective theory. This leads 
to well known effects such as the time dilatation correction [41], but can also account for interference contributions 
with processes where no top pairs are produced, but which have the same final state [44] (see Fig. 6). The gauge 
cancellation already mentioned for class (a) also applies here and leads to the cancellation of interactions among the 
top quarks and their decay products caused by ultrasoft gluons for the total cross section at NLL order [45, 46] and 
even at NNLL order [44] . In Ref. [44] it was also shown that within the nonrelativistic effective theory the imaginary 
matching conditions can lead to ultraviolet phase space divergences that require additional renormalization. A 
complete treatment of all the class (c) electroweak effects at NLL has not yet been achieved. 

A crucial prediction of the Higgs mechanism is that the Higgs Yukawa coupling to quarks X q is related to the quark 
masses by m q = X q V. At the e + e~ Linear Collider the top quark Yukawa coupling can be measured from top quark 
pair production associated with a Higgs boson, e + e~ — > tiH. This process is particularly suited for a light Higgs 
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Figure 6: (a) Cuts in Standard Model diagrams contributing to absorptive parts in the matching conditions for the tt effective 
theory currents. Within the effective theory the absorptive parts describe the interference of the diagrams shown in (b). 
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Figure 7: Typical constellation of momenta for the process e + e — > ttH in the large Higgs energy endpoint region. 



boson since the cross section can then reach the 1-2 fb level. Assuming an experimental precision at the percent level, 
QCD and electroweak radiative corrections need to be accounted for in the theoretical predictions. The Born cross 
section was already determined some time ago in Ref. [47]. For the 0(a s ) QCD one-loop corrections a number of 
references in various approximations exist [48-50]. On the other hand, the full set of one-loop electroweak corrections 
was obtained in Refs. [51, 52] and also in Ref. [53]. In Ref. [52] a detailed analysis of various differential distributions 
of the cross section er(e+e~ — ► ttH) can be found. For high energies these fixed order predictions are sufficient to 
reach the required theoretical precision. 

In Ref. [54] the phase space region where the Higgs energy is large and the tt pair becomes collinear to balance the 
large Higgs momentum was analyzed (see Fig. 7). It was found that since in this kinematic region the ti invariant 
mass approaches 2m t , it is governed by nonrelativistic dynamics in analogy to the physics relevant for the threshold. 
Thus, the usual fixed-order treatment breaks down and the nonrelativistic effective theory description known from 
the ti threshold has to be applied. For high cm. energies above about 700 GeV the large Higgs energy endpoint 
region is very small, and a usual fixed-order treatment is sufficient for the theoretical determination of the total 
cross section. However, for smaller cm. energies (as available at the first phase of the ILC running), or for larger 
Higgs masses, the endpoint region increases with respect to the full available phase space and an effective theory 
treatment for the endpoint region becomes mandatory. In Ref. [54] a QCD endpoint analysis in the framework of 
the effective theory vNRQCD [55] was carried out. The results sum singular terms in the endpoint region cx (a s /v) n 
and (a s lnv) n at NLL order, where v is the cm. velocity of the top quarks in the ti cm. frame. Typical results for 
the Higgs energy spectrum for small cm. energies showing the fixed-order QCD predictions and the LL and NLL 
order nonrelativistic effective theory predictions in the large Higgs energy region are shown in Fig. 8. The vertical 
lines mark the Higgs energy where the top velocity v = 0.2, which roughly divides the regions where fixed-order and 
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Figure 8: Higgs energy spectrum at LL (dotted lines) and NLL (solid lines) order in vNRQCD for the soft renormalization 
scales us = (0.1, 0.2, 0.4)m t and at the Born level and at 0(a s ) for fi = -^/s, for the parameters (a) yfs — 500 GeV, 

m H = 120 GeV and (b) y/s = 550 GeV, m H = 140 GeV. The top IS mass has been set to m t is = 180 GeV and the other 
parameters are T t = 1.55 GeV, M z = 91.1876 GeV, M w = 80.423 GeV, c w = M w /M z and a" 1 = 137.034. 



v/£ [GeV] 


m H [GeV] 


cr(Born) [fb] 


a(a 3 ) [fb] 


cr(NLL) [fb] 


cr(NLL) 
<r(Bom) 


er(NLL) 

<r(tx„) 


500 


120 


0.151 


0.263 


0.357(20) 


2.362 


1.359 


550 


120 


0.984 


1.251 


1.342(37) 


1.364 


1.073 


550 


160 


0.134 


0.207 


0.254(12) 


1.902 


1.226 


600 


120 


1.691 


1.939 


2.005(30) 


1.185 


1.034 


600 


160 


0.565 


0.700 


0.745(18) 


1.319 


1.065 


700 


120 


2.348 


2.454 


2.485(13) 


1.058 


1.012 


700 


160 


1.210 


1.303 


1.328(11) 


1.098 


1.020 



Table III: Cross sections and K factors for o"tot(e e 
mf = 180 GeV. 



ttH) for various cm. energies and Higgs masses and top quark mass 



effective theory computations are valid. In Table III the impact of the summation of the Coulomb singularities and 
the logarithms of the top quark velocity in the endpoint region is analyzed numerically for various choices of the 
cm. energy and the Higgs mass. For all cases the top quark mass ml s = 180 GeV is used and the other parameters 
are fixed as in Fig. 8. Here, er(Born) refers to the Born cross section and a(a s ) to the 0(a s ) cross section in fixed- 
order perturbation theory using fi = y^s as the renormalization scale. The term er(NLL) refers to the sum of the 
0(a s ) fixed-order cross section for v > 0.2 using fi = yfsv and the NLL nonrelativistic cross section for \v\ < 0.2 with 
the soft vNRQCD renormalization parameter /ig = 0.2wt. The summation of endpoint singularities is particularly 
important for smaller cm. energies, since it leads to an additional significant enhancement of the total cross section 
in a region where the cross section is rather small. This enhancement could improve the prospects of top- Yukawa 
coupling measurements at the first phase of the ILC (see section 2.5). 



2.5. Towards a Precise Measurement of the Top Quark Yukawa Coupling at the ILC [56] 

A. Juste 

The top quark Yukawa coupling (At) is the largest coupling of the Higgs boson to fermions. A precise measure- 
ment of it is very important since it may help unravel the secrets of the Electroweak Symmetry Breaking (EWSB) 
mechanism, in which the top quark could possibly play a key role. For rrih < 2mt, a direct measurement of At is 
possible via associated tth production, both at the LHC and a future ILC. At the LHC, the expected accuracy [57] 
is £)At/A t ~ 12 — 15% for rrih ~ 120 — 200 GeV, assuming an integrated luminosity of 300 fb -1 . Existing feasibility 
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studies at the ILC [58] predict an accuracy of <5A t /A t <~ 6 — 10% for rrih ~ 120 — 190 GeV, assuming y/s = 800 GeV 
and 1000 fb -1 . However, currently the baseline design for the ILC only contemplates a maximum center-of-mass 
energy of 500 GeV. Therefore, it is very relevant to explore the prospects of this measurement during the first phase 
of the ILC, especially in view of the limited accuracy expected at the LHC; for a number of years, the combination 
of results from the LHC and ILC will yield the most precise determination of A t . 

A preliminary feasibility study at y/s = 500 GeV was performed in Ref. [59], which we briefly overview here. 
Indeed, the measurement of At at y/s = 500 GeV is more challenging than at y/s — 800 GeV. On the one hand, 
the reduced phase-space leads to a large reduction in a t ih (e.g. °"^ r " — 0.16(2.5) fb at y/s = 500(800) GeV, for 
rrih = 120 GeV). On the other hand, the cross section for many background processes is significantly increased. This 
analysis assumed rrih = 120 GeV and focused on the tih — ► (£vb)(jjb)(bb) decay channel (BR ~ 30%). The dominant 
background is tijj, followed by tiZ, although other non- interfering backgrounds (e.g. LF+VF - ) were also considered. 
Signal and backgrounds were processed through a fast detector simulation. After basic preselection cuts, the signal 
efficiency was found to be ~ 50% and the S : B ~ 1 : 100. In order to increase the sensitivity, a multivariate analysis 
using a Neural Network (NN) with 23 variables was performed. The final selection consisted of an optimized cut 
on the NN distribution. Assuming an integrated luminosity of 1000 fb -1 , the expected total numbers of signal and 
background events were 11 and 51, respectively, resulting in (5\t / \t) stat — 33%. Based on previous experience [58], 
the addition of the fully hadronic decay channel was expected to ultimately lead to {6X t /X t ) stat ~ 23%. While this 
analysis is already rather sophisticated, significant improvements are expected from the usage of a more efficient 
6-tagging algorithm or a more optimal treatment of the kinematic information. In the next sections we discuss 
additional sources of improvement which are currently under investigation. 

The precise measurement of At requires accurate theoretical predictions for o~tih- Currently, one- loop QCD and 
electroweak corrections are available. However, at y/s = 500 GeV and for rrih > 120 GeV, the kinematic region 
where the ti system is non-relativistic dominates. As discussed in Ref. [54] and Section 2.4, in this regime Coulomb 
singularities are important and need to be resummed within the framework of the vNRQCD effective theory, leading 
to large enhancements factors in the cross section relative to the Born level. At the ILC, because of ISR and 
beamstrahlung (BS), the event-by-event center-of-mass energy (y/1) will be lower than the nominal one, thus bringing 
the ti system deeper into the non-relativistic regime. In order to compute the expected a t t h including these effects, 
the 11-fold Born differential cross section for e + e~ — > tih — > W + bW~bh was multiplied by a K-factor defined as 
K(Eh,Vs) = (daS^ L /dEh)/(daS? rn /dEh), where Eh stands for the Higgs boson energy in the e + e~ rest-frame, 
and then folded with ISR and BS structure functions. The NLL differential cross section was kindly provided by 
the authors of Ref. [54]. Fig. 9 (left and center) compares the Born (for off-shell top quarks) and NLL differential 
cross sections for different values of \/§, assuming m\ s = 180 GeV and rrih = 120 GeV. The ratio of these two 
curves defines K(Eh,V§) and can be significantly larger than 1, especially for low values of Since the NLL 
prediction is only valid for Eh < E^ ax (where E"^ ax effectively corresponds to a cut on the top quark velocity in the 
ti rest-frame of j3 t < 0.2), we currently set K(Eh, V§) = 1 for Eh > E h nax , although in practice, it should be possible 
to use K(E h ,\Ts) = {da° lh as) / dE h ) / (da^° h rn / dE h ). Table IV compares the predicted Born and "NLL-improved" 
a t ih for different scenarios, illustrating the large impact of radiative effects in the initial state. This underscores the 
importance of being able to predict these effects to the percent level. While the impact of ISR cannot be reduced, it 
might be possible to find an optimal operating point of the accelerator, as far as this measurement is concerned, in 
terms of BS and total integrated luminosity. Finally, it is found that, for rrih = 120 GeV, resummation effects can 
increase u t i h by a factor of ~ 2.4 with respect to the Born cross section used in the previous feasibility study. 

So far, all feasibility studies of this measurement have assumed unpolarized beams. Currently, the baseline design 
for the ILC only includes longitudinal electron beam polarization (\P e - \ — 0.8). Positron beam polarization (\P e + | ~ 
0.6) is considered as an option. The ratio of the polarized cross section (for arbitrary longitudinal beam polarization) 
(o~p _p + ) to the unpolarized cross section (er ) is given by up _p /<tq = (1 — P e -P e +)(1 — P e ffALR), where 
P e ff = (P e - — P e +)/(l — P e -P e +) denotes the "effective polarization" and Alr is the "left-right asymmetry" of the 
process of interest [60]. Therefore, two potential enhancement factors can in principle be exploited: the first one 
requires having both beams polarized, the second one requires A^r / and a judicious choice of the signs of P e - and 
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Table IV: Comparison of the Born and NLL a t t h for different scenarios regarding radiative effects in the initial state. 



(ISR,BS) 


&tth (fb) (Born) 


°tth (fb) ("NLL-improved") 


Enhancement factor 


(off,off) 


0.157(f) 


0.357(2) 


2.27 


(off,on) 


0.f06(f) 


0.252(3) 


2.38 


(on, on) 


0.0735(8) 


0.179(2) 


2.44 




E h (GeV) E b (GeV) 



Figure 9: Left and center: comparison of the Born (dashed) and NLL (solid) do t i h / dEh for different values of i/s, assuming 
m] s = 180 GeV and ran = 120 GeV. The dotted line indicates the value of E™ ax . Right: ratio of polarized to unpolarized 
cross section for different values of (P e -,P e +). 

P e + in order to have P^^Alr < 0. In the case of SM tih production, Alr ~ 0.44, essentially independent of y/s in the 
range ~ 0.5 — 1.0 TeV. Assuming (Alr)sm, Fig. 9 (right) shows the cross section enhancement factor as a function 
of P e +, for different values of P e - . The optimal (realistic) operating point would be (P e -,P e +) = (—0.8, +0.6), 
achieving an increase in o~ t ih by a factor of ~ 2.1 with respect to the unpolarized case. Unfortunately, this choice 
does not help reduce the dominant background, which is increased by a similar factor. Nevertheless, the net result is 
still an improvement in the statistical precision on A t by ~ 45%, which would be an argument in favor of including 
positron polarization in the baseline design. For (P e - , P e +) = ( — 0.8, 0), only a modest increase in a t ih by a factor of 
~ 1.3 would be achieved. It is important to realize that, in order to choose the sign of P e - and P e + , it is necessary 
to know the sign of Alr. Anomalous couplings in the ttj and ttZ vertices could possibly lead to deviations in Alr 
from the SM prediction. Unfortunately, due to the limited precision in the measurement of the ttZ couplings [12, 16], 
the LHC is not expected to provide any useful constraints on the sign of Alr. Therefore, at the ILC the first step 
should be to perform measurements of the polarized er t j in order to determine the sign of Alr, and thus fix the 
signs of P e - and P e + (the magnitudes should be the largest possible). On the other hand, the measurement of At 
requires a percent-level and model-independent determination of the tl/y and tiZ couplings, which typically benefits 
from changing the beam polarization. Therefore, it would be desirable to optimize the running strategy to maintain 
the largest possible cr t i hl needed for a precise measurement of A t , while meeting the precision goals for measurements 
of top quark couplings. 

We have studied the prospects of a precise measurement of the top quark Yukawa coupling during the first phase 
of the ILC. Taking into consideration an existing feasibility study, and the additional enhancement factors to 
discussed here, we anticipate a precision of (6Xt/Xt)stat ~ 10% for mu = 120 GeV, assuming y/s = 500 GeV and 
1000 fb" 1 . 
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2.6. The tt Threshold at an e+e~ Collider 

Y. Kiyo 

It is well known that QCD corrections to tt production near threshold develop a Coulomb singularity. The structure 
of the cross section is 



o~tot(e + e — > tf) = a Be 



i , - (<±\ . c ( 2 ) ( + c (3) f £!• y + . . . + c («) f «• 



(7) 



where a Born is the Born cross section and v = \/l — 4m 2 / s is the velocity of the top quarks. In Eq. 7, we have 
pulled out the Coulomb singularity (a s /v) n explicitly so that the nth coefficient starts with O(v ), 

c {n) = c („,0) + c (n,l) v + c (n,2) y 2 + . . . _ (g) 

Near threshold the kinematics of the top quarks is non-relativistic and v <~ a s holds. To get a meaningful cross 
section we thus have to sum up the Coulomb singularities <~ a s /v ~ 1. The leading order (LO) cross section then 
contains the Coulomb singularities <~ (a s /v) n to all orders of the coupling expansion, 



LO 



„ Ec K0) M\ 

n=0 ^ ' 



^„VcWp . (9) 



The c™' 1 terms are suppressed by v compared with the LO and result in next-to-leading order (NLO) corrections. 
These terms are known since long. Next-to-next-to leading order (NNLO) calculations have been completed some 
time ago by several groups [61]. They sum up all the corrections of c'"' 2 ) v 2 (a s /v) n to all orders n. 




Figure 10: Typical QCD loop diagram which yields the Coulomb singularity ~ (a s /v) n , n — the number of gluon exchanges. 

The Coulomb singularity originates from potential gluons, which have a typical momentum fc M ~ mt(v 2 ,v), ex- 
changed between almost on-shell top and anti-top quarks in Fig. 10. For the kinematics of potential gluons, the 
corresponding propagator reduces to the Coulomb potential 

% = (10) 

q 2 

where CV = 4/3. In the gluon propagator, the energy component of q was set to zero because of the potential gluon's 
kinematics. Using this gluon propagator in place of normal gluon propagators in the figure, one may reproduce the 
c O,o) coe fg c i en t s to all orders in Eq. 8. This procedure can be systematically extended to higher orders using effective 
field theory (EFT) techniques. The constructed effective field theories are versions of non-relativistic QCD called 
pNRQCD/vNRQCD [55]. We do not go into details of the EFT but summarize several features. 

• The tt pair at threshold is created by the production current J = Cj [tat\ + ■ ■ ■ in the EFT, where Cj is the 
Wilson coefficient of the corresponding current in the EFT, and the dots denote subleading operators. 

• The produced top quark pair forms bound state resonances by exchanging potential gluons; the binding Coulomb 
potential is Vq = VoVc(q), where Vo = 1 + 0{a s ) includes QCD corrections to the leading order Coulomb 
potential from loop diagrams. 

• There are subdominant potentials, e.g. V\ = Vig 2 /|q|, which have to be taken into account if one wants to go 
beyond the NLO cross section. 



PLEN0043 



2005 ALCPG & ILC Workshops - U.S.A. 



• A yet new type of corrections is known at next-to-next-to-next-to leading order (NNNLO), which is referred to 
as ultra-soft corrections because of the typical gluon momentum of k° ~ |fc| ~ mv 2 . 

An important point here is that the EFT scheme is a systematic way to sum up higher order corrections, and it 
makes practical calculations easier compared to those of full QCD. Currently complete NNLO total cross sections 
are known in both (semi-) analytical and numerical ways. The results are summarized in Ref. [61]. 

In the following we briefly survey recent attempts in going beyond NNLO and in resumming potentially large 
logarithms for the threshold cross section in the EFT framework. 

In the EFT calculation of the tt cross section, higher order QCD corrections enter through the higher order 
coefficients Vi of potentials and through subdominant potentials, e.g. Vi(q). Furthermore, there are corrections from 
ultra-soft gluons, starting to contribute at NNNLO. They cannot be written in form of potentials, as they include a 
noninstantaneous, dynamical propagation in time, while the potentials are all instantaneous. The first correction to 
Vo is referred to as a Coulomb correction, 

v (») = x + ^(0oL(q) + f) + •••+ (^)" • (group factors, L(q), etc.), (11) 

where /3q = 11 — (2/3)ny is the coefficient of the QCD (3— function, L(q) — ln(/i 2 /q 2 ), and /i is the QCD renor- 
malization scale. The QCD corrections to all the coefficients can be found in Ref. [62] up to 0(a^), except the 

~ (3") (3) ~ 

third order constant term, so-called CI3 in the literature. Using Vq — Vq Vq, the Coulomb correction to the total 
cross section was computed in Ref. [63]. The other corrections to the tt cross section at NNNLO due to subleading 
potentials and ultra-soft gluons are not known yet. In Fig. 11 we show the tt threshold total cross section including 
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Figure 11: Normalized cross section R = a(e + e~ —> tt)/a(e + e~ — > fi + n~) in the threshold regime including the resummation 
of Coulomb corrections. Figure from Ref. [63] . 

the Coulomb corrections as calculated in Ref. [63]. The curves are obtained using analytical results of successive 
Coulomb corrections from LO to NNNLO, and the line denoted by "NNNLO exact" is obtained by numerically 
solving the Schrodinger equation using Vq at NNNLO. The difference between NNNLO and "NNNLO exact" is 
yet of higher order as the numerical method contains rigid NNNLO plus additional contributions which come from 
iterations of the potential Vq. Studying the renormalization scale dependence of these predictions, the theoretical 
uncertainty for the NNNLO Coulomb corrections was estimated to be about 5% of the cross section. Although the 
full corrections at NNNLO are not known, the IS energy of the tt resonance, which corresponds to the peak position 
of the cross section, is known analytically at NNNLO from Refs. [62, 64]. The result shows that the uncertainty in 
the choice of the scale is small and the perturbative series for the IS energy is stabilized (see Ref. [65]). 

In higher order calculations there appear potentially large logarithms of ratios of scales. In the threshold regime 
we have large logarithms of v in the coefficients c^" ,J \ arising from ratios of the largely different scales of the problem, 
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namely the energy of the top quarks, E t ~ mv 2 , the momentum p ~ m t v of almost on-shell top quarks, and the 
cm. energy yfs ~ 2mt- The origin of these logarithms is related to UV divergences in the EFT. They can be 
resummed using renormalization group (RG) arguments [55, 66], resulting in RG improved cross section predictions, 
o~ll, o~nll, etc. The RG improvement was extensively studied recently, and the next-next-leading-log (NNLL) cross 
section was calculated in Ref. [67] (apart from the NNLL running of a current for which the anomalous dimension 
is only fully known at NLL, see also Ref. [68, 69]). The RG improved cross section then contains all terms of order 
Enm( a s/")"( a s' ln ') m - The resummation of the logarithms leads to a reduction of the scale dependence of the 
normalization from 20% to about 6% (at NNLL), especially at energies around the IS 1 peak. This is illustrated in 
Fig. 12, where the scale dependence of the fixed order (left panel) and RG improved (right panel) cross sections are 
shown for different orders. 

The different studies discussed above strengthen our confidence that the estimate of the uncertainty of the top 
quark mass determination from Ref. [61], Am t < 100 MeV, is realistic and will not be spoiled by uncalculated higher 
order corrections. There is even hope to believe that, if different approaches turn out to converge, the estimated error 
may shrink. For determinations of the top quark total width, the top Yukawa coupling y t and the strong coupling 
constant from the threshold scan theoretical uncertainties are comparable to the expected experimental errors. An 
ultimate theoretical normalization error of at most 3% is desirable. 




O.o t — , — i — , — i — , — i — , — i — , — i — , — i — , — i — , — j o.O t — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — j 

346 347 348 349 350 351 352 353 354 346 347 348 349 350 351 352 353 354 

yfs(GeV) yfs(GeV) 



Figure 12: Scale dependence of the fixed order cross section (a) and the RG improved cross section predictions (b) in different 
orders. Each band of curves is obtained by varying the vNRQCD renormalization scale, choosing v — 0.15,0.2,0.3. Figures 
from Ref. [68]. 



2.7. Differential Distributions and Experimental Aspects of the tt Threshold Scan 

T. Teubner 

The most precise determination of the top quark mass will come from a dedicated threshold scan at s/s ~ 2mt- 
This is mainly a consequence of the fact that by performing a counting experiment of color singlet tt pairs one 
can avoid most of the systematic uncertainties inherent in a kinematic reconstruction of the decay products of the 
colored quarks. These uncertainties, which are closely related to the problem of the mass definition, will ultimately 
limit the determination of nit as done at the Tevatron and soon at the LHC. Detailed studies have shown that via 
a threshold scan, measurements with very small statistical and systematic errors will be possible at the ILC [2]. 
A multi-parameter fit for the top quark mass rnt, width T t , the strong coupling a s and the top Yukawa coupling 
resulted in experimental errors of about Amj ~ 20 MeV, AT t ~ 30 MeV, Aa s ~ 0.0012, depending on details of 
the fit. (At threshold, the top Yukawa coupling can only be measured for a light Higgs and even then with less than 
30% accuracy.) Such a precision will only be achieved if the accuracy of the theoretical predictions can match the 
experimental one, and recently a lot of effort has been invested in further improvements of the theory (see sections 
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2.4 and 2.6). Although most of the information in the fit will come from the precise measurement of the total tt cross 
section, differential distributions are needed for several reasons. 

• Experimentally, cuts are needed to discriminate the signal form backgrounds, so the measured cross section 
can never be the fully inclusive total cross section. 

• Distributions are required to build realistic (higher order) Monte Carlo generators for the signal process. 

• Using additional observables beyond otot adds information, helping to disentangle correlations among the 
parameters determined from a threshold scan, and increases the sensitivity to possible New Physics in top 
production and decay. 

Differential distributions used so far are the top quark momentum distribution da/dpt and the forward-backward 
asymmetry Afb of the cross section, but theoretical studies exist also for the polarization of the top quarks. In the 
following we will briefly discuss these observables and their role in the threshold scan, and comment also on the 
issue of rescattering corrections. 

0.04 
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Figure 13: Top quark momentum distribution in the threshold regime. 

The top quark momentum distribution do~/dp t is available in the framework of the next-to-next-to leading order 
(NNLO) calculations [70] (see also Ref. [71]). Close to threshold its form resembles a 15 Coulomb-like wave function, 
becoming more symmetric well above threshold. Fig. 13 shows do~/dp t at a fixed cm. energy of 349 GeV for two 
different masses m t = 174.5, 175.5 GeV, where the width of the band is obtained by varying a s between 0.115 and 
0.121. The peak position of da/dpt is not much dependent on a s but is very sensitive to mt, and a change of 
a s mainly affects the normalization of the distribution. This is in contrast to 0t o t, wnere m t and cx s are strongly 
correlated parameters when fitting theoretical predictions to (simulated) data. 

Interference of the leading vector current (through 7 and Z exchange) with the suppressed axial vector contribution 
(from Z only) leads to a forward-backward asymmetry Afb ■ The size of the asymmetry depends on how much the 
corresponding 5 and P wave resonance contributions overlap and hence on the top quark width Tt, but less on a s 
(and mi), see Fig. 14. 

As already demonstrated in experimental studies [2], the two observables da/dp t and Ape can be used together 
with a tot to make best use of the available information and to disentangle the correlations between the parameters 
mt,T t ,a s (,yt) in a multi-paramter fit. However, the latest theoretical developments such as NNNLO corrections, 
renormalization group improvement by summing large logarithms at NNLL order (see section 2.6) or the effect of 
EW corrections (see section 2.4) are dealing with the total cross section, whereas distributions are available at fixed 
NNLO only. This situation is quite common, as higher order corrections are increasingly difficult for differential 
cross sections, and not many distributions are known beyond NLO. Nevertheless, from a pragmatic point of view 
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Figure 14: Dependence of the forward-backward asymmetry Afb on the top quark width (left) and the strong coupling a s 
(right panel). 



and for use in Monte Carlos, it is legitimate to use the best available prediction for the total cross section together 
with distributions available at lower order. By rescaling the distributions using the total cross section, their scale 
dependence will decrease and a consistent normalization will be ensured. 




Figure 15: Feynman diagram responsible for rescattering corrections between the top and the b from the t decay. 



All recent NNLO, NNNLO as well as the renormalization group improved calculations are not taking into account 
the so-called "rescattering corrections" , i.e. interactions between the t, t quarks and the b, b quarks from the top 
decays, or between b and 6 (see Fig. 16 for a typical Feynman diagram). However, such corrections were calculated 
numerically in Refs. [72, 73] to NLO accuracy. They are strongly suppressed (vanishing at NLO) in the inclusive cr to t 
but generally important for distributions and typically of order 10%. They therefore should be included in a realistic 
Monte Carlo description. 

In addition to differential distributions, also the top quark polarization has been studied in the threshold 
regime [72-74] . Predictions are available at NLO accuracy for all three polarization components, including rescatter- 
ing corrections. Even for unpolarized e + , e~ beams the top quarks are (longitudinally) polarized to 40%, and for the 
foreseen beam polarization the top quarks will be highly polarized. Due to the short life-time the top polarization will 
be transmitted undisturbed to the decay b's and W's and allow very interesting studies. Making use of the polariza- 
tion, observables can be defined which are sensitive to the top's electric dipole moment (probing CP-violation beyond 
the Standard Model) or to anomalous top couplings like V + A admixtures to the weak interaction. However, more 
realistic experimental studies are needed in order to fully explore the potential of polarization in the threshold regime. 
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Experimental aspects of the threshold scan 

Precision measurements from a ti threshold scan require a very precise knowledge of both the average cm. energy 
{\/s} and of the luminosity spectrum dL/d^/s [75]. This is not an easy task, as the beam dynamics at a linear collider 
is not as constrained as it was e.g. at LEP, and only a single measurement of bunches before collision will be possible. 
The average cm. energy can be determined by energy (up- and downstream) spectrometers in dedicated beam line 
inserts, discussed in Working Group 4 of the Snowmass Workshop. A possible problem here could be a bias between 
the spectrometer measurements and the collision cm. energy. Another possibility to measure the beam energy is 
the use of physics processes like Z pair production, or radiative return (through photon radiation) calibrated to the 
Z peak [76]. 

The luminosity spectrum is determined by (a) the beam spread, (b) beamstrahlung, and (c) initial state radiation 
(ISR). All three effects will lead to a smearing of the ti threshold cross section, resulting in a significant reduction of 
the effective luminosity and hence the observed cross section, 

CT ° bs (Vi) = -^ ( L{x)a{xyfs)dx. (12) 

L Q Jo 

The influence of the three effects is demonstrated in Fig. 16. The beam spread will typically be ~ 0.1% at the 
ILC and will cause comparably little smearing (though additional beam diagnostics may be required to measure and 
monitor the beam spread), but beamstrahlung and ISR are very important. The luminosity spectrum dL/d^/s will 
lead to a systematic shift in the extracted top mass which must be well understood; otherwise it could become the 
dominant systematic error. In the past, ti threshold studies were carried out under the assumption that the spectrum 
is basically known. However it has turned out that the precise determination of dL/d^/s is a challenging task. The 
proposed method is to analyse the acollinearity of (large angle) Bhabha scattering events, which is sensitive to a 
momentum mismatch between the beams but insensitive to the absolute energy scale [77]. For this, the envisioned 
high resolution of the forward tracker will be very important to achieve the required accuracy. 3 




vJi [GeV] 



Figure 16: Left: Smearing of the theoretical ti cross section ('default') by beam effects and initial state radiation. Right panel: 
Simulation of beam spread, beamstrahlung and ISR as distributions of x — \fs/ \/so (where ^/so is the nominal cm. energy 
of the machine). 

In the simulation of the luminosity spectrum beamstrahlung is the main unknown. Integrated ILC accelerator 
simulations for these machine dependent effects are becoming available, including effects from the linac beam spread, 
the beam diagnostics and delivery system and the simulation of the collision dynamics using the package Guinea-pig. 
At this workshop Stewart Boogert presented first results of spectrum parameterizations (using the package CIRCE), 



2 This section can only provide a brief review and the reader interested in more details is referred to the presentations of Stewart 
Boogert at this workshop. 

3 An additional complication is, that in Bhabha scattering ISR is not completely factorizable from final state radiation. 
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based on new simulations from G. White for ^s = 350 GeV. The new parameterizations show variations from 
earlier results which were based on beam simulations at higher energies only. The parameterizations will be used 
as input in the simulation of the Bhabha scattering (based on the package BHWIDE) for the acollinearity analysis. 
These involved simulations are required to extract dL/dyfs and to understand the uncertainties of the reconstructed 
spectrum. Only then can one quantify in physics analyses like the tt threshold scan, to which extent the uncertainties 
in the luminosity spectrum affect the accuracy of the intended measurements like the top quark mass, width, etc. 

The effects due to ISR are in principle machine independent and calculable to high precision in QED. Nevertheless, 
ISR effects complicate the measurement of beamstrahlung and accelerator energy spread, and one has to ensure that 
the precision of the theoretical formulae used in the Monte Carlo codes is sufficient. Similarly, the accuracy of 
differential distributions from programs used for the simulation of the wide angle Bhabha events must be assessed, 
and recent theoretical results [78] may have to be included. 

The recent developments discussed here make it clear that much work remains to be done before the sophisticated 
analysis techniques for the tt threshold scan are understood at the required level of accuracy. For a complete study, 
which should take into account the total cross section, distributions, full beam effects including asymmetric boosts 
and detector effects, a full Monte Carlo generator will be needed, and a project for this has been started at Snowmass. 
With a more detailed experimental study it will also be possible to optimize the scan strategy, i.e. how the available 
luminosity should be distributed among different energy scan points. 

The analysis of beamstrahlung and its impact on the physics program at different energies will also be relevant for 
the optimization of the machine design. Concerning detectors, the demands of the top pair counting experiment are 
most probably not problematic for the existing detector designs. However, the requirements for Bhabha scattering and 
radiative Z return might set the requirements of the low angle tracking system and the electromagnetic calorimeter. 
Finally, it should also be noted, that the understanding of other physics processes like the WW or SUSY thresholds 
will largely benefit from the top case, and tt at threshold should be regarded as the benchmark. 

3. QCD EFFECTS IN ILC PHYSICS 

3.1. Precision QCD at the ILC: e+e~ — > 3 Jets 

A. Gehrmann-De Ridder, T. Gehrmann, E.W.N. Glover 

The production of light quark-antiquark pairs in electron-positron annihilation gives rise to final states containing 
QCD jets. Depending on the amount of additional hard QCD radiation, one obtains final states with a certain 
number of jets: if only quark and antiquark are hard, two-jet final states are produced, one additional hard gluon 
yields a three-jet final state, two extra gluons or a secondary quark-antiquark pair can give rise to four-jet final 
states and so on. Studying these multi-jet final states, one can probe many aspects of perturbative QCD. Three-jet 
final states and related event shape observables were studied extensively at LEP in order to determine the strong 
coupling constant a s , which controls the probability of radiating a hard gluon in these events. The measurement of 
four-jet-type observables at LEP established the gauge group structure of QCD. Five-jet and higher multiplicities 
were sometimes considered in new physics searches, where QCD-induced processes form a theoretically calculable 
background. 

The LEP measurements of three-jet observables are of a very high statistical precision. The extraction of a s from 
these data sets relies on a comparison of the data with theoretical predictions. Comparing the different sources of 
error in this extraction [79], one finds that the purely experimental error is negligible compared to the theoretical 
uncertainty. There are two sources of theoretical uncertainty: the theoretical description of the parton-to-hadron 
transition (hadronization uncertainty) and the theoretical calculation of parton-level jet production (perturbative 
or scale uncertainty). Although the precise size of the hadronization uncertainty is debatable and perhaps often 
underestimated, it is certainly appropriate to consider the scale uncertainty as the dominant source of theoretical 
error on the precise determination of a s from three-jet observables. This scale uncertainty arises from truncating the 
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perturbative QCD expansion of jet observables at the next- to-leading order (NLO) and can be improved considerably 
by computing next-to-next-to-leading order (NNLO) corrections. 

Given the planned luminosity of the ILC, one expects that this collider will deliver jet-production data of a statis- 
tical quality similar to LEP. An attractive perspective of such a measurement at the ILC would be to determine the 
evolution of a s over a wide energy range, which is potentially sensitive to new physics thresholds. Concerning uncer- 
tainties on such a determination, it is worthwhile to note that the hadronization corrections become less important 
at higher energies, thus leaving the scale uncertainty as dominant source of theoretical error. Experimental aspects 
of such measurements, especially issues related to the beam energy profile (which were irrelevant at LEP) were not 
studied up to now, and certainly deserve further attention. 

In the recent past, many steps towards the NNLO calculation of e + e~ — ► 3 jets have been accomplished (see 
Ref. [80] and references therein). Foremost, the relevant two- loop 1 — > 3 matrix elements are now available. One- loop 
corrections to 1 — > 4 matrix elements have been known for longer and form part of NLO calculations of e + e~ — » 4 jets. 
These NLO matrix elements naturally contribute to e + e~ — ► 3 jets at NNLO if one of the partons involved becomes 
unresolved (soft or collincar). In this case, the infrared singular parts of the matrix elements need to be extracted 
and integrated over the phase space appropriate to the unresolved configuration to make the infrared pole structure 
explicit. As a final ingredient, the tree level 1 — > 5 processes also contribute to e + e~ — > 3 jets at NNLO. These contain 
double real radiation singularities corresponding to two partons becoming simultaneously soft and/or collinear. To 
compute the contributions from single unresolved radiation at one-loop and double real radiation at tree-level, one 
has to find subtraction terms which coincide with the full matrix elements in the unresolved limits and are still 
sufficiently simple to be integrated analytically in order to cancel their infrared pole structure with the two-loop 
virtual contributions. In the following, we present a new method, named antenna subtraction, to carry out NNLO 
calculations of jet observables and discuss its application to e + e~ — ► 3 jets. 

In electron-positron annihilation, an m-jet cross section at NLO is obtained by summing contributions from (m+1)- 
parton tree level and m-parton one-loop processes: 



do- N LO = / (&°~NLO _ da NLo) + \ 



d(7 NLO + 



d( ?NLO 



The cross section da^ LO is the (m+ l)-parton tree-level cross section, while dcrJ^ LO is the one- loop virtual correction 
to the m-parton Born cross section dcr B . Both contain infrared singularities, which are explicit poles in 1/e in da^ LO , 
while becoming explicit in da^ LO only after integration over the phase space. In general, this integration involves 
the (often iterative) definition of the jet observable, such that an analytic integration is not feasible (and also not 
appropriate). Instead, one would like to have a flexible method that can be easily adapted to different jet observables 
or jet definitions. Therefore, the infrared singularities of the real radiation contributions should be extracted using 
infrared subtraction terms. One introduces do~f[ LO , which is a counter-term for da^ LO , having the same unintegrated 
singular behavior as da^ LO in all appropriate limits. Their difference is free of divergences and can be integrated 
over the (m + l)-parton phase space numerically. The subtraction term daff LO has to be integrated analytically over 
all singular regions of the (m + l)-parton phase space. The resulting cross section added to the virtual contribution 
yields an infrared finite result. Several methods for constructing NLO subtraction terms systematically were proposed 
in the literature [81-84]. For some of these methods, extension to NNLO was discussed [86] and partly worked out. 
We focus on the antenna subtraction method [81, 82], which we extend to NNLO. 

The basic idea of the antenna subtraction approach at NLO is to construct the subtraction term derjy i0 from 
antenna functions. Each antenna function encapsulates all singular limits due to the emission of one unresolved 
parton between two color-connected hard partons (tree- level three-parton antenna function). This construction 
exploits the universal factorization of phase space and squared matrix elements in all unresolved limits, depicted in 
Fig. 17. The individual antenna functions are obtained by normalizing three-parton tree-level matrix elements to the 
corresponding two-parton tree-level matrix elements. 

At NNLO, the m-jet production is induced by final states containing up to (m + 2) partons, including the one-loop 
virtual corrections to (m + l)-parton final states. As at NLO, one has to introduce subtraction terms for the (m+1)- 
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Figure 17: Illustration of NLO antenna factorization representing the factorization of both the squared matrix elements and 
the (m + l)-particle phase space. The term in square brackets represents both the antenna function and the antenna phase 
space. 




Figure 18: Illustration of NNLO antenna factorization representing the factorization of both the squared matrix elements and 
the (m + 2)-particle phase space when the unresolved particles are color connected. The term in square brackets represents 
both the antenna function and the antenna phase space. 



and (to + 2)-parton contributions. Schematically the NNLO m-jet cross section reads, 

do- NNLO = / (dcr NNLO - da NNLO ) + / dcr NNLi 

+ / [do- NNLO - dtT NNLo) + / d(T NNL 

Jd<S> m+1 v ' Jd* m+ i 



d$„ 



NNLO 

' Jd* m+ i 

, V,2 

ao ~NNLO ' 



where da NNLO denotes the real radiation subtraction term coinciding with the (TO + 2)-parton tree level cross section 
d^NNLO m a ^ singular limits [87]. Likewise, da N ^ l LO is the one-loop virtual subtraction term coinciding with the 
one-loop (to + l)-parton cross section da^ NLO in all singular limits [88]. Finally, the two- loop correction to the 
TO-parton cross section is denoted by da N ' NLQ . 

Both types of NNLO subtraction terms can be constructed from antenna functions. In da NNLO , we have to distin- 
guish four different types of unresolved configurations: (a) One unresolved parton but the experimental observable 
selects only m jets; (b) Two color-connected unresolved partons (color-connected); (c) Two unresolved partons that 
are not color connected but share a common radiator (almost color- unconnected); (d) Two unresolved partons that 
are well separated from each other in the color chain (color-unconnected). Among those, configuration (a) is prop- 
erly accounted for by a single tree-level three-parton antenna function like used already at NLO. Configuration (b) 
requires a tree-level four-parton antenna function (two unresolved partons emitted between a pair of hard partons) 
as shown in Fig. 18, while (c) and (d) are accounted for by products of two tree-level three-parton antenna functions. 

In single unresolved limits, the one- loop cross section da N ^ LO is described by the sum of two terms [88]: a 
tree-level splitting function times a one-loop cross section and a one-loop splitting function times a tree-level cross 

VS 1 

section. Consequently, the one-loop single unresolved subtraction term da NI j LO is constructed from tree-level and 
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Figure 19: Illustration of NNLO antenna factorization representing the factorization of both the one-loop "squared" matrix 
elements (represented by the white blob) and the (m+l)-particle phase space when the unresolved particles are color connected. 

V S 1 

one- loop three-parton antenna functions, as sketched in Fig. 19. Several other terms in da N] ^ LO cancel with the 
results from the integration of terms in the double real radiation subtraction term daff NLO over the phase space 
appropriate to one of the unresolved partons, thus ensuring the cancellation of all explicit infrared poles in the 
difference do%„ LO - da^ S / LQ . 

Finally, all remaining terms in do~fj NLO and da^^LO have to be integrated over the four-parton and three-parton 
antenna phase spaces. After integration, the infrared poles are rendered explicit and cancel with the infrared pole 
terms in the two- loop squared matrix clement da N ' NL0 . 

The subtraction terms da^ LO , daf !NLO and dcr]^£ require three different types of antenna functions corre- 
sponding to the different pairs of hard partons forming the antenna: quark-antiquark, quark-gluon and gluon-gluon 
antenna functions. In the past [81, 82], NLO antenna functions were constructed by imposing definite properties 
in all single unresolved limits (two collinear limits and one soft limit for each antenna). This procedure turns out 
to be impractical at NNLO, where each antenna function must have definite behaviors in a large number of single 
and double unresolved limits. Instead, we derive these antenna functions in a systematic manner from physical ma- 
trix elements known to possess the correct limits. The quark-antiquark antenna functions can be obtained directly 
from the e+e~ — > 2j real radiation corrections at NLO and NNLO [89]. For quark-gluon and gluon-gluon antenna 
functions, effective Lagrangians are used to obtain tree-level processes yielding a quark-gluon or gluon-gluon final 
state. The antenna functions are then obtained from the real radiation corrections to these processes. Quark-gluon 
antenna functions were derived [90] from the purely QCD (i.e. non-supersymmetric) NLO and NNLO corrections to 
the decay of a heavy neutralino into a massless gluino plus partons [91], while gluon-gluon antenna functions [92] 
result from the QCD corrections to Higgs boson decay into partons [93] . 

All tree-level three-parton and four-parton antenna functions and three-parton one-loop antenna functions are 
listed in Ref. [80], where we also provide their integrated forms, obtained using the phase space integration techniques 
described in Ref. [94]. 

In Refs. [80, 95] we derived the l/7V 2 -contribution to the NNLO corrections to e+e~ — ► 3 jets. This color factor 
receives contributions from 7* — > qqggg and 7* — > qqqqg at tree-level [96], 7* — > qqgg and 7* — > qqqq at one-loop [97] 
and 7* — ► qqg at two- loops [98]. The four-parton and five-parton final states contain infrared singularities, which 
need to be extracted using the antenna subtraction formalism. 

In this contribution, all gluons are effectively photon-like, and couple only to the quarks, but not to each other. 
Consequently, only quark-antiquark antenna functions appear in the construction of the subtraction terms. 

Starting from the program EERAD2 [81], which computes the four-jet production at NLO, we implemented the 
NNLO antenna subtraction method for the 1/N 2 color factor contribution to e + e~ — > 3j. EERAD2 already contains 
the five-parton and four-parton matrix elements relevant here, as well as the NLO-type subtraction terms. 

The implementation contains three channels, classified by their partonic multiplicity: (a) in the five-parton channel, 
we integrate da^ NLO — daff NLO ; (b) in the four-parton channel, we integrate dcrJ^ LO — dcr]^^ ; (c) in the three- 
parton channel, we integrate da^^LO + ^nnlo + ^^nlo ■ ^ ne numerical integration over these channels is carried 
out by Monte Carlo methods. 
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By construction, the integrands in the four-parton and three-parton channel are free of explicit infrared poles. 
In the five-parton and four-parton channel, we tested the proper implementation of the subtraction by generating 
trajectories of phase space points approaching a given single or double unresolved limit. Along these trajectories, 
we observe that the antenna subtraction terms converge locally towards the physical matrix elements, and that the 
cancellations among individual contributions to the subtraction terms take place as expected. Moreover, we checked 
the correctness of the subtraction by introducing a lower cut (slicing parameter) on the phase space variables, and 
observing that our results are independent of this cut (provided it is chosen small enough). This behavior indicates 
that the subtraction terms ensure that the contribution of potentially singular regions of the final state phase space 
does not contribute to the numerical integrals, but is accounted for analytically. 

As a final point, we noted in Ref. [80] that the infrared poles of the two- loop (including one- loop times one- loop) 
correction to 7* — > qqg are canceled in all color factors by a combination of integrated three-parton and four-parton 
antenna functions. This highly non-trivial cancellation clearly illustrates that the antenna functions derived here 
correctly approximate QCD matrix elements in all infrared singular limits at NNLO. They also outline the structure 
of infrared cancellations in e+e~ — > 3j at NNLO, and indicate the structure of the subtraction terms in all color 
factors. 

In this talk, we discussed the theoretical prerequisites for performing precision QCD studies on existing LEP data 
and at the ILC. In particular, the precise extraction of the strong coupling constant a s requires improved theoretical 
predictions to reduce the scale error inherent to calculations in perturbative QCD. At present, this extraction relies 
on the calculation of e+e~ — ► 3 jets at NLO accuracy, and we reported on progress towards the NNLO calculation. 

This calculation requires a new method for the subtraction of infrared singularities which we call antenna sub- 
traction. We introduced subtraction terms for double real radiation at tree level and single real radiation at one 
loop based on antenna functions. These antenna functions describe the color-ordered radiation of unresolved partons 
between a pair of hard (radiator) partons. All antenna functions at NLO and NNLO can be derived systematically 
from physical matrix elements. 

Using this method, we implemented the NNLO corrections to the subleading color contribution to e+e~ — ► 3 jets 
into a flexible parton level event generator program, and are currently proceeding with the implementation [99] of 
the full set of color factors relevant to the NNLO corrections to e+e~ — > 3 jets. 

3.2. Study of V L V L -> ti at the ILC Including 0(a 2 a ) Corrections [100] 

S. Godfrey 

Understanding the mechanism of electroweak symmetry breaking (EWSB) is a primary goal of the LHC and ILC 
[57]. While much effort has been devoted to the weakly interacting weak sector scenario the strongly interacting weak 
sector (SIWS) remains a possibility. Because the i-quark mass is the same order of magnitude as the scale of EWSB 
it has long been suspected that t-quark properties may provide hints about the nature of EWSB and the subprocess 
VlVl — * ti has been suggested as a probe. While VlVl — * ti can be studied at both hadron colliders and e + e~ 
colliders, the overwhelming QCD backgrounds will likely make it impossible to study the VlVl — » ti subprocess at 
the LHC [101]. In contrast, the ILC offers a much cleaner environment. But to be able to attach meaning to precision 
measurements it is necessary to understand radiative corrections, both electroweak and QCD. Here we show 0{a s ) 
corrections to the tree level electroweak VlVl — > ti process in the SM at the ILC. Due to space limitations we point 
the interested reader to Ref. [102] for a more detailed account and a more complete set of references. 

We are interested in the subprocesses VV — ► ti which occur in the processes e + e~ — > £1^2 + VV — > £1^2 + ti where 
£1(2 is vv for the W + W~ — ► ti subprocess and e + e~ for the ZZ — ► ti subprocess. The vector bosons are treated 
as partons inside the e + and e~ using the effective boson approximation [103, 104]. The total cross section is then 
obtained by integrating the W (or Z) luminosities with the subprocess cross section [105]. 

The 0(a s ) corrections for the processes W + W~ —* ti ZZ — > ti are calculated using the FeynArts, FormCalc and 
LoopTools packages [106]. The QCD corrections to W + W~ — > ti are shown in Fig. 20 (left side). The infrared 
singularity in the vertex corrections are canceled by the soft contributions from the process W + W~ — ► tig which 
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Figure 21: (left side) The K-factors as a function of yfs e+e - for e + e~ — ► vvtt (via W^W^ fusion). The solid line is 
for Mh = 120 GeV, the dashed line for Mh = 500 GeV, the dotted line for Mh = 1 TeV, and the dot-dashed line for 
Mh = oo (LET), (right side) The K-factor as a function of Mh for e + e~ — > vvtt (via W^W^ fusion). The solid line is for 
\fs e+e - — 500 GeV, the dashed line for y/s e + e - = 1 TeV, the dotted line for \fs e + e - = 2 TeV, and the dot-dashed line for 
\fs e+e - = 3 TeV. See text for an explanation of the K-factor. 



are shown in Fig. 20 (right side). We regulate the IR-singularity by introducing a gluon mass which is equivalent 
to standard dimensional regularization for processes with no triple gluon vertex present. This approach has the 
additional benefit that varying the value of the gluon mass acts as a check of the numerical cancellations between 
the different contributions. 

We include in our results the kinematic cuts m t i > 400 GeV and > 10 GeV. Since the longitudinal scattering 
cross section is much larger than the TT and TL cases and it is the longitudinal gauge boson processes which 
corresponds to the Goldstone bosons of the theory we will henceforth only include results for VlVl scattering. 

The QCD corrections to longitudinal scattering are often presented as a K-factor, normally defined as the ratio 
of the NLO to LO cross sections. Because the 0(a s ) QCD corrections we calculated are LO corrections to a tree 
level electroweak result we take the K-factor to be the ratio of the cross section with the 0{a s ) QCD corrections and 
the tree level electroweak cross sections. The K-factors for a(e + e~ — ► vvtt) which goes via W^W^ fusion and for 
cr(e + e~ — > e + e~tt) is shown in Fig. 21 (left side) as a function of y/s e + e -- The 0(a s ) QCD corrections are largest 
for M H = 500 GeV with K-factors ranging from over 1.2 for \/i e+e - = 500 GeV to 1.15 for \/i e+e - = 1 TeV. The 
corrections decrease as y/s e+e - increases. The variation of the K-factor with Mh is shown in Fig. 21 (right side). 
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The fact that the K- factor is largest for M H = 500 GeV in Fig. 21 (left side) and that it peaks at M H ~ 400 GeV 
in Fig. 21 (right side) is a threshold effect which is an artifact of the kinematic cut we imposed on the tt invariant 
mass. The important point is that the QCD corrections are not insignificant compared to the effects we might wish 
to study such as top Yukawa couplings or anomalous VVti couplings. 

3.3. QCD at a Photon Collider [107] 

Z. Sullivan 

A terascale photon collider will provide a unique opportunity to understand the resolved hadronic structure of light 
beyond a few GeV. This structure completely dominates the QCD cross section if the invariant mass of hadronic final 
state particles is less than ~ Wmax/3. This is clear from Fig. 22(a), where we see the cross section for bb production 
as a function of invariant mass M^. The first calculation of the uncertainty in this cross section is presented in 
detail in the Proceedings of this workshop [107]. The result, the dashed band surrounding the upper solid curve in 
Fig. 22(b), is that the cross section cannot be predicted to better than a factor of 5. There is no way to improve 
this prediction without measuring the gluon, charm, and bottom parton distributions for the photon in situ at this 
collider. 

Loosening the cuts used for Fig. 22(a) will provide an extremely clean sample of events with a gluon in the initial 
state. Hence, the gluon structure should be quickly well-measured. The peak near 90 GeV is from a Z resonance, 
for which 40% of the events come from a cc initial state. Therefore, a clean extraction of the charm structure can be 
made by reconstructing the Z peak in multiple channels (particularly Z — > /i + /i~). Finally, the long resolved-resolved 
tail above the Z peak is almost entirely due to bb + bb + bb collisions. Therefore, given enough data, even the bottom 
structure may be accessible. 

The study [107] focused on a photon-photon collider, but another option being considered is a photon-electron 
collider. In general, it is more difficult to cleanly extract the gluon PDF in a 7— e collision than a 7-7 collision, 
because the cross section at high invariant mass (> 20 GeV) is smaller, the decay products tend to be boosted more 
forward into less-well instrumented regions of the detector, and an additional deconvolution must be performed to 
remove the effect of extracting an almost-real photon from the electron. Nevertheless, this option should be examined 
in more detail as it may be simpler to construct a 7-e collider. 
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4. NEW PHYSICS EFFECTS ON TOP QUARK PROPERTIES 



4.1. Top Compositeness at Colliders 



K. Agashe 



Consider the Randall- Sundrum (RSI) model [108] which is a compact slice of AdSs 



ds 2 = e-^^-q^dx^dXv+rlde 2 , - 



ir<9<ir 



(13) 



where k is the curvature scale and the extra-dimensional interval is realized as an orbifolded circle of radius r c . The 
two orbifold fixed points, 9 = 0, ir, correspond to the "UV" (or "Planck) and "IR" (or "TeV") branes respectively. 
In warped spacetimes the relationship between 5D mass scales and 4D mass scales (in an effective 4D description) 
depends on location in the extra dimension through the warp factor, e~ fc ' 6 '' rc . This allows large 4D mass hierarchies 
to naturally arise without large hierarchies in the defining 5D theory, whose mass parameters are taken to be of order 
the observed Planck scale, Mpi <~ 10 18 GeV. For example, the 4D massless graviton mode is localized near the UV 
brane while the Higgs sector is taken to be localized on the IR brane. In the 4D effective theory one then finds 



A modestly large radius, i.e., knr c ~ log (Mp;/TeV) ~ 30, can then accommodate a TeV-size weak scale. Kaluza- 
Klcin (KK) graviton resonances have masses m <~ fce~ fc7rr <\ These masses are at the TeV-scale, since their wave 
functions are also localized near the IR brane. 

In the original RSI model, it was assumed that the entire SM was localized on the TeV brane, and that only 
gravity propagated in the full 5D space. Thus, the effective UV cut-off for gauge and fermion fields and hence the 
scale suppressing higher-dimensional operators is at a TeV, the same scale which sets the Higgs sector. However, 
bounds from electroweak precision tests (EWPT) on this cut-off are approximately 5-10 TeV, while those from flavor 
changing neutral currents (FCNCs) such as K — K mixing are around 1000 TeV. Stabilizing the electroweak scale 
thus requires fine-tuning; even though RSI explains the big hierarchy between the Planck and electroweak scales, it 
has a "little" hierarchy problem between the weak scale and the TeV scale cut-off. 

An attractive solution to this problem is to allow the SM gauge [109] and fermion [110, 111] fields to propagate in 
the extra dimensional bulk. We first explain how bulk fcrmions enable us to evade flavor constraints. The localization 
of the wavefunction of the massless chiral mode is controlled by the 5d mass term for each fermion, which in units of 
k is denoted by the c-parameter. In the warped scenario, for c > 1/2 (c < 1/2) the zero mode is localized near the 
Planck (TeV) brane, whereas for c = 1/2, the wave function is flat. We therefore choose c > 1/2 for light fcrmions so 
that the effective UV cut-off at the location of the light fermions is much greater than a TeV, suppressing dangerous 
FCNCs. This naturally results in a small 4D Yukawa coupling to the Higgs on the TeV brane without any hierarchies 
in the fundamental 51? Yukawa couplings [110-112]. Similarly, we choose c <C 1/2 for the top quark to obtain an 
0(1) Yukawa coupling. We can also show that in this scenario with bulk gauge fields high-scale unification of gauge 
couplings can be accommodated. 

Since gauge fields are also in the bulk, their excited KK modes induce additional effects in flavor physics and 
in EWPT, which are calculable in the 5D effective field theory. For example, the couplings of KK modes to light 
fermions are flavor-dependent, giving FCNCs. However, this flavor dependence is small. KK modes, just like the 
Higgs, are localized near the TeV brane, whereas the light fermions are near the Planck brane. The FCNCs are 
therefore proportional to the Yukawa couplings, resulting in a suppression of flavor violation. This scenario has an 
analog of the Glashow-Iliopoulos-Maiani (GIM) mechanism of the SM, resulting in the suppression of the calculable 



Early studies showed that the corrections to EWPT from gauge KK modes are too large, unless the KK mass is 
greater than 10 TeV [113]. Such a large KK scale results in a little hierarchy problem between the weak and KK 
scales. The localization of the light fcrmions near the Planck brane reduces the contribution of gauge KK modes to 



Weak Scale ~ M Pi e 



-k-!rr c 



(14) 



FCNCs [111, 112]. 
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Figure 23: Shift in the coupling of tn to Z as a function of KK mass in GeV on horizontal axis and cr on vertical axis for the 
choice cl = 0.4. The solid, long-dashed and short-dashed lines correspond to shifts of 5%, 10% and 15%, respectively. 



two observables - the S parameter and 4-fermion operators, but the observable called the T parameter still gives 
stringent constraints. 

In Ref. [114], it was shown that this problem can be avoided by extending the electroweak gauge group in the 
bulk to SU{2)l x SU(2)r x U(1)b-l- Such an extension provides a custodial isospin symmetry to protect the T 
parameter from large corrections. Thus, KK masses as low as 3 TeV are allowed by oblique EW data, significantly 
ameliorating the little hierarchy problem. 

We now consider the top and bottom quarks. It is clear that we prefer c « 1/2 for ti, to obtain a top Yukawa 
coupling of 0(1) without too large a 5D Yukawa coupling, but this implies a large shift in coupling of to the 
Z-boson unless the KK scale is larger than a TeV. This shift occurs due to the large coupling of b^ to the KK Z 
modes, which mix with the zero-mode Z via the Higgs vev [114]: 



5 



9*i 



m| 



l-2c fl / kwr c 5~2c R \ 



+ 7To ^ • (15) 



KK gauge (OAItukk) 2 3 - 2c R V 2 '*(3-2c R )J 

Here, itikk ~ 2.45 ke~ km ' c is defined to be the mass of the lightest gauge KK mode. Thus, there is a tension 
between obtaining the top Yukawa and not shifting the coupling of b^ to Z . As a compromise, KK masses ~ 5 TeV 
are consistent with a shift in g b z L of approximately 0.25% for cl <~ 0.4. The tR must therefore be localized near the 
TeV brane: cr ~ 0. Such a profile for tR leads to sizable shift in its coupling to the Z via exchange of KK Z-bosons, 
and also via KK ti, modes. The additional contribution from KK exchange is due to zero- mode tR mixing with 
KK th via a Higgs vev which then couples to the Z, giving the shift 

6 (ffz ) „ 1/2 / mt ^l/2-c L \ 2 



9z R 



KK fcrmion ^— ' -2/3 sin 9 W \ m (») 



L 



Here m t a/1/2 — ci is the mass term coupling zero- mode and KK top quarks, and the KK th masses are given by 
m. ( n) w 7T ke~ k7Tr -(n - c L /2) w 0.78 m KK (n - c L /2). 

The total shift in the coupling of tR is plotted in Fig. 23. As seen in the figure, we obtain an 10% shift for KK 
masses of a few TeV. Smaller KK masses are not allowed by EWPT, and larger masses lead to large fine-tuning. 
Observability of an effect of this size might be difficult at the LHC since the sensitivity of the LHC is only at the 
20% level for a shift in the axial coupling of the top quark to the Z. It should be possible at the ILC, which has 
sensitivity at the few percent level for shifts in both axial and vector couplings of the top quark to Z [11]. There are 
similar shifts in couplings of Higgs to W/Z due to its profile being localized near TeV brane just like for tR. 
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Figure 24: The corrections to the ttZ axial and vector couplings in the original Littlest Higgs model (left panel) and its 
extension with T-parity (right panel). The regions in which the ILC would observe no deviation from the SM are shaded. 



Finally, an intriguing aspect is that via the AdS/CFT correspondence [115], such a scenario is conjectured to be 
dual to a purely AD theory with a composite Higgs boson [116], with the light fermions being elementary and tn 
being composite. This provides an intuitive understanding for the large shifts in the couplings of tn and Higgs to Z. 
Hence, these signals might be valid for general composite Higgs models as well. 



4.2. Top quark properties in Little Higgs Models [118] 

C.F. Berger, M. Perelstein, F. Petriello 

In this section, we study the corrections to the top quark properties in "Little Higgs" models of electroweak 
symmetry breaking [117], and compare the expected deviations from the SM predictions with expected sensitivities 
of experiments at the LHC and the ILC. In the Little Higgs models, electroweak symmetry is driven by the radiative 
effects from the top sector, including the SM-like top and its heavy counterpart, a TeV-scale "heavy top" T. Probing 
this structure experimentally is quite difficult. While the LHC should be able to discover the T quark, its potential 
for studying its couplings is limited [119, 120]. Direct production of the T will likely be beyond the kinematic reach 
of the ILC. However, we will show below that the corrections to the gauge couplings of the SM top, induced by 
its mixing with the T, will be observable at the ILC throughout the parameter range consistent with naturalness. 
Measuring these corrections will provide a unique window on the top sector of the Little Higgs. 

Little Higgs models contain a light Higgs boson which is a composite of more fundamental degrees of freedom. A 
generic composite Higgs model must become strongly coupled at an energy scale around 1 TeV, leading to unaccept- 
ably large corrections to precision electroweak observables. In contrast, Little Higgs models remain perturbative until 
a higher energy scale, around 10 TeV. The hierarchy between the Higgs mass and the strong coupling scale is natural 
and stable with respect to radiative corrections. Because of the special symmetry structure of the theory, the Higgs 
mass vanishes at tree level, as do one-loop quadratically divergent diagrams. The mass term is dominated by the 
logarithmically divergent one-loop contribution from the top quark, which triggers electroweak symmetry breaking. 

Many Little Higgs models have been proposed in the literature. We will consider two examples in this study, the 
"Littlest Higgs" model [121], and its variation incorporating T parity [122]. We will study the effects on the tiZ 
vertex in these models; for a more detailed study, see Ref. [118]. 

Corrections to the gauge couplings of the top quark in Little Higgs model arise from two sources: the mixing 
of the (left-handed) top with the heavy top T, and the mixing of the SM gauge bosons W ± , Z° with their heavy 
counterparts, and W#. Using the superscripts "t" and "g" to denote the contributions from these two sources, 
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the corrections to the ttZ coupling can be written as 

\ 2 T v 2 g\ 



Sg z R e = 



r 2 H 2 



4/ 



2 r.1 



S 9l S = Jp 



2.9 



A s t/j 



9 Z n- 



r 2 s 2 



(17) 



Here, gf R are the SM left- and right-handed ttZ couplings, gy = (g R + crf)/2 and g\ = (gf, — gf)/2 are their 
vector and axial combinations, cw, sw are respectively the cosine and sine of the weak mixing angle, and = sin?/;, 
a,/} = cos-0- In the original Littlest Higgs model [121], both the gauge sector shift and the top sector shift occur; 
in the T-parity model [122], only the top sector shift is present. The predicted shifts in the ttZ axial and vector 
couplings for mr = 0.5,1.0, and 2.0 TeV, and At = 0.5,1,2, are plotted in Fig. 24 (left panel), along with the 
experimental sensitivities expected at the LHC [12] and the ILC [11]. The mixing angle ip is varied between and 
7r/2. Note that the shifts have a definite sign. While only a rather small part of the parameter space is accessible at 
the LHC even with 3000 fb _1 integrated luminosity, the ILC experiments will be able to easily observe the shifts in 
most of the parameter space preferred by naturalness considerations (however, the prospects for observation at the 
LHC improve when additional final states such as bb + 4 j are included; see Ref. [124]). Similarly, shifts in the Wtb 
coupling can be probed via deviations in the top quark width at the ILC [118]. 



4.3. Testing CPT Symmetry with Top Quark Physics [125] 

J.A.R. Cernbranos 

The viability to observe evidence of CPT violation in the top sector has been analyzed through the measurement 
of a mass difference between top and anti-top [125]. This study has been focused on the CPT violating ratio of the 
top quark, RcPT(t) = 2(m* — mi)/(mt + Tn^). The present constraints from the Tevatron are approximately 10%, 
and they could be reduced by one order of magnitude at the LHC or ILC. The most promising studied channel is the 
lepton plus jets channel for top anti-top production. However, other techniques to reconstruct the top mass could 
also be very interesting, such as the analysis of the J/psi from b decay at the LHC, which improves the systematic 
uncertainties. Single top production could also be studied, since a combination of different measurements would be 
necessary in order to consider CPT violation as the explanation of any exotic data. 

Di-lcpton channel: Di-lepton events originating predominantly from tt — > W + (-^ l + v) bW~(-^ l~v) b, with £ — e 
or fj,, have been used in Tevatron to measure the top quark pole mass supposing an identical mass for the top and 
anti-top quarks. The same data can be used to study CPT violation through a double peak in the reconstructed 
invariant mass associated to the lepton and b quark coming from the single decay of the top or anti-top [125]. By 




Figure 25: Schematics of the top and anti-top decays in the dilepton channel. 



using the Tevatron data accumulated at Fermilab from 1992 through 1995 [126], it is possible to find the bound 
\Rcpr(t)\ < 0.13 at the 95% C.L. [125]. On the other hand, the sensitivity of the LHC can be estimated as 
Rc'PT(t) = 0.03 at the 95% c.l. following an analogous analysis [125]. 

Lepton plus jets channel: A more promising signal is provided when one of the W decays leptonically while the 
other one decays hadronically: tt — > W + {-^ £ + vb)bW~ qq')b. In fact, the inclusive lepton plus jets channel 
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provides a larger and cleaner sample of top quarks, whose mass can be reconstructed directly using the hadronic part 
of the decay. The invariant mass of the three jets coming from the top (rrijjb = fnj q jgj b ) or anti-top (rrijjb = m j q j 9 j i ) 
presents a peak at the top (mt) or anti-top (mt) mass respectively. 

The estimate combining the CDF [127] and D0 data [128] gives a more constraining bound of Rcpr(t) < 9.2 x 10~ 2 . 
The sensitivity of the LHC is also better in this channel since both statistical and systematic uncertainties are expected 
to be improved. Indeed, the LHC will be able to test the CPT violation of the top quark to almost one order of 
magnitude better than the present constraints: \RcPT(t)\ — 0.014 at the 95% C.L. or equivalently, m t — mi ~ 2.4 
GeV [125]. 




Figure 26: Schematic example of the top and anti-top decays in the lepton plus jets channel. 

The same analyses can be performed with the ILC, where an increase of the statistical uncertainties but a decrease 
of the systematic ones is expected [129]. The importance of these last uncertainties leads to a small improvement of 
the sensitivity in relation to the LHC. 
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